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ABSTRACT 

We use a N-body/hydrodynamic simulation to forecast the future encounter be- 
tween the Milky Way and the Andromeda galaxies, given current observational con- 
straints on their relative distance, relative velocity, and masses. Allowing for a com- 
parable amount of diffuse mass to fill the volume of the Local Group, we find that 
the two galaxies are likely to collide in a few billion years - within the Sun's lifetime. 
During the the interaction, there is a chance that the Sun will be pulled away from 
its present orbital radius and reside in an extended tidal tail. The likelihood for this 
outcome increases as the merger progresses, and there is a remote possibility that 
our Sun will be more tightly bound to Andromeda than to the Milky Way before the 
final merger. Eventually, after the merger has completed, the Sun is most likely to be 
scattered to the outer halo and reside at much larger radii (> 30 kpc). The density 
profiles of the stars, gas and dark matter in the merger product resemble those of 
elliptical galaxies. Our Local Group model therefore provides a prototype progenitor 
of late-forming elliptical galaxies. 

Key words: galaxy:evolution — galaxies:evolution — galaxies: formation — galax- 
iesdnteractions — Local Group — methods :N-body simulations. 



1 INTRODUCTION 

It is well known that the Milky Way (MW) and Andromeda 
(M31) are the two largest members of the Local Group of 
galaxies. Together with their ~ 40 smaller companions, the 
Milky Way and Andromeda comprise our galactic neigh- 
borhood, and as such, represent the nearest laboratory, and 
therefore the most powerful tool, to study the formation and 
, evolution of galactic structure. 

Like most extragalactic groups, the Local Group is 
very likely to be decoupled from the cosmological expan- 
sion and is now a gravitationally bound collection of galax- 
ies. This notion is supported by the observed relative mo- 
tion between its two largest galaxies; namely, the Milky 
Way and A ndromeda are moving tow ard each other at ~ 
120 km s~ 1 (|Binnev fc Tremainelll987h . Unfortunately, this 
motion alone does not indicate whether the Local Group 
is bound or not. The unknown magnitude of Andromeda's 
transverse velocity adds uncertainty into the present day or- 
bital parameters and therefore the past and future evolution 
of the Local Group. 

Barring the uncertain transverse velocity of An- 
dromeda, a considerable amount of information can be in- 
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ferred about the Local Group pro vided a plausible set o f 
assumptions. Nearly 50 years ago iKahn fc Woltierl (|l959h 
pioneered the "timing argument," in which the Milky Way 
and Andromeda are assumed to form within close proximity 
to each other, during the dense early stages of the Universe, 
before they were pulled apart by the general cosmological 
expansion. They have subsequently reversed their path and 
are approaching one another owing to their mutual gravi- 
tational attraction. According to the timing argument, the 
Milky Way and Andromeda have now traced out nearly a full 
period of their orbital motion which is governed by Kepler's 
laws. By assuming that the system has no angular momen- 
tum, and given the current separation, velocity of approach, 
and the age of the Universe, the timing argument yields es- 
timates for the mass of the Local Group (> 3 x 10 12 M Q ), 
the semi-major axis of the orbit (< 580 kpc), and the 
time of the next close pa ssage (> 4 Gyr) (see Sec. 10.2 of 
iBinnev fc Tremaind[l987l) . 

While the seminal results of IKahn fc Woltierl (|l959h 
were an early indication of the large mass-to-light ratio in 
the Local Group and therefore the presence of dark matter, 
they also began a nearly five decade long quest to under- 
stand the past, present, and future of our Local Group. In 
particular, a number of studies have extended the original 
timing argument by allowing for various angular momenta, 
by including more realistic or time-dependent mass distri- 
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butions, by adding the effects of mass at scales beyond that 
of the Local Group, or testing its validity using numerical 
simu l ations (see, e.g ., Peeble s et al.l 1989; iFich fc Tremaine! 
199ll;IValtonen et al II 19931; iPeeblesll 1994 [Peebles et al.l | 2001 



Sawa fc Fuiimotoll2005l: lLoeb et al.ll2005l ; iLi fc WhiteH2007 
van der Marel fe Guhathakurtall2007l ). 



One of the most intriguing developments stemming 
from the various studies of the Local Group is an estimate of 
the transverse velocity of Andromeda. By employing the ac- 
tion principle to the motions o f galaxies with i n and near 
(< 20 Mpc) the Local Group, IPeebles et ail l|200ll ) con- 
cluded that the transverse velocity of Andromeda is less 
than 200 km s " 1 . Using the well mea sured transverse ve- 
locity of M33 l|Brunthaler et all [2005) and numerical sim- 
ulations that tracked the potential tidal disruption during 
M33's past encounters with Andromeda. lLoeb et alJ (|2005 ) 
found an even smaller estimate, ~ 100 km s _1 , for the trans- 
verse velocity. While future astrometric observations using 
SIM0 and GAIA0 will be able to accurately measure the 
proper motion of Andromeda, the low values favored by 
these papers suggests that the Local Group is indeed a grav- 
ftationally bound system. 

Provided that the Local Group is gravitationally bound, 
and that the Milky Way and Andromeda are heading to- 
wards each other, one must admit the possibility that they 
will eventually interact and merge. This outcome appears 
inevitable given the massive halos of dark matter that likely 
surround the Milky Way and Andromeda. Numerical exper- 
iments have robustly concluded that dark matter halos can 
exert significant dynamical friction, and are sponges that 
soak up energy and angular momentum leading to a rapid 
merger l|Barneslll988h . 

Even though the eventual merger between the Milky 
Way and Andromeda is common lore in Astronomy, the 
merger process has not been addressed by a comprehen- 
sive numerical study . The one exception is a paper by 
iDubinski et al.l l|l996l ) that presented a viable model for 
the Local Group and numerically simulated the eventual 
merger between the M ilky Way and Andromeda. However, 
IDubinski etafl (|l996l ) utilized this Local Group model and 
its numerical evolution to study the production of tidal tails 
during such an encounter and the possibility to use the struc- 
ture of this tidal materi al to probe the dark matter poten- 
tial. While the study by IDubinski et all (|l996l ) provided the 
first enticing picture of the future encounter between the 
Milky Way and Andromeda, (for a more recent and hig her 
resolution version of this simulation, see lDubinskill2006h . it 
was neither designed to detail the merger dynamics includ- 
ing intergalactic material, nor outline the possible outcomes 
for the dynamics of our Sun, nor quantify properties of the 
merger remnant. In addition, the last decade has produced 
a number of improved models for the structure of the Milky 
Way and Andromeda as well as the properties of the intra- 
group medium. 

In this paper we quantitatively predict when the inter- 
action and merger of the Milky Way and Andromeda will 
likely occur and forecast the probable dynamics of the Sun 
during this event. We achieve this goal by constructing a 



model for the Local Group in Sj2]that satisfies all observa- 
tional constraints. We then evolve this model using a self- 
consistent N-body/hydrodynamic simulation, as described 
in f|3] The generic properties of the merger, including the 
merger timescale, the possible evolution of our Solar Sys- 
tem, and properties of the merger remnant, are outlined in 
Sj4] Finally, we conclude in SJU 



2 A MODEL OF THE LOCAL GROUP 

The distribution of mass within our Local Group of galaxies 
has been a long-standing question in astrophysics. It is clear 
that much of the matter is associated with the two largest 
galaxies in the Local Group: the Milky Way and Andromeda. 
Moreover, these two spiral galaxies are likely to be embedded 
in an ambient medium of dark matter and gas. 



2.1 The Milky Way and Andromeda 

There are a number of different models for both the Milky 



l nere are a number ot airrerent models tor potn tne Milky 
Way and Andromeda galaxies (see, e.g. , Klypin et al.l |2002| ; 
IWidrow fc Dubinskill2005l ; ISeigar et al.ll2006l . and reference 
therein). These studies generally enlist a myriad of obser- 
vational data to infer the distribution of baryons, while the 
dark matter, which dominates the gravitational potential, 
is set to match distrib utions extracted from cosmological N- 
body simulations fe.g. jNavarro et al.lll99r]| ). Together, these 
models specify the total mass distribution out to the virial 
radius (~ 200 - 300 kpc). 

In our model of the Local Group we start by adopting 
the models for the Milky Way and Andromeda favored by 
iKlypin et all (|2002t ). Within these models, the baryons are 
contained entirely within the rotationally supported expo- 
nential disk and central bulge. These components are then 
surrounded by a massive dark-matter halo, which has nearly 
20 times the mass as the baryons, as specified by the mass 
fractions, mj and m^, defined as the bulge and disk mass, re- 
spectively, divided by the total mass. The exponential disk, 
of radial disk scale radius Rd, also contains a set fraction / 
of its mass in collisional gas that can cool and form stars. 
Both the bulge and d ark h alo components are assumed to 
follow the iHernquistl (|l990h profile. The bulge scale radius 
a is fixed to be 20% of the radial disk scale radius Rd. The 
dark-matter profile is defined by its concentration c, spin 
parameter A, and total virial mass M200 and virial circu- 
lar velocity V200 (at the radius r2oo where the average inte- 
rior density is 200 times the critical cosmic density today, 
rhocrit = 10 -29 g cm -3 ), which are all listed in Table [1] 
The numerical construction of these models employs meth- 
ods com monly used to construct equilibrium disk galaxies 



see e.g..lHernquisdll993al;ISpringel fc Whitell 19991 ; ISpringell 



1 http:/ /planetquest. jpl.nasa.gov/SIM/ 

2 http://www.sci.esa.int/gaia/ 



(sec, e.g., ricrnqmst tyyja; springel & Whi 
2000;; ICox et al.ll200rj ; ISpringel et al.l l2005) 



2.2 The Orbit 

Given the adopted parameters of the two largest galaxies 
in the Local Group, we must now define their orbital pa- 
rameters and any ambient medium in which the system will 
be embedded. There are a few empirical constraints that 
must be considered. First, at the current epoch, the sepa- 
ration between the Milky Way and Andromeda is 780 kpc 
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Table 1. Properties of the Milky Way (MW) and Andromeda 
(M31) models used in the work (see text for definitions). 



Property 


MW 


Andromeda 


V200 (km s" 1 ) 


145 


170 


M200 (10 12 M Q ) 


1.0 


1.6 


c 


12 


12 


\ 

A 


U.UOl 


U.UoO 


™d 


0.041 


0.044 




2.2 


3.6 


f 


0.3 


0.3 


m b 


0.008 


0.012 


a 


0.4 


0.7 


N dm 


475,500 


755,200 




14,350 


24,640 


N 9<ls 


6,150 


10,560 




4,000 


9,600 



ijMcConnachie et all I200I iRibas et all l2005t ). Second, the 
Milky Way and Andromeda are approaching each other at a 
radial speed of 120 km s~ , ass uming a local circular spee d 
of 220 kms _1 (see Sec. 10.2 of iBinnev fc Tremaind 1 19871 ). 
These observational facts tightly constrain any dynamical 
model of the Local Group since their fractional error bars 
are estimated to be less than ~ 5%. 

Less well constrained is the current estimate 
for the transverse velocity of Andromeda. As 



tioned in JT] sp eeds of < 200 
recent models dPeebles et al 



Ivan der Marel fc Guhathakurta 



men- 

km s -1 are favo r ed by 
12051 ; lLoeb et al.l 120051 ; 
20071 ). but depend upon 
assumptions regarding the distribution of mass within the 
Local Group and its initial state. We will therefore gauge 
the success of our Local Group model by its ability to re- 
produce the above three observations, but it should be kept 
in mind that the first two observations have significantly 
less leeway than the third. 

Constraints on the spin orientation of both the Milky 
Way and Andromeda with respect to the orbital plane of 
the merger originate from the current position and orienta- 
tion of Andromeda in the night sky. While these details are 
necessary for a complete model of the Local Group, they do 
not significantly influence the timing of the merger between 
the Milky Way and Andromeda. The spin orientation will, 
however, affect the disk morphology during the merger, a 
detail that will be addressed when we attempt to track the 
possible dynamical fate of our Sun in H4.3.2I 

2.3 The Local Intragroup Medium 

One plausible starting model for the Local Group is to 
follow the logic originally e mployed by the timing argu- 
ment (|Kahn fc Woltierlll959r ). i.e., that the mass of the Lo- 
cal Group is entirely contained within the Milky Way and 
Andromeda and their motion is a simple two body prob- 
lem governed by Kepler's equations. In practice, however, 
most of t he recent implementations of the timing argument 
(see, e.g..lBinnev fc Tremainelll987l ; iFich fc Tremaing|l99ll ; 
iLi fc Whitel2007l ) generally yield masses for the Local Group 
(> 3 x 10 12 Mq) that exceed the total masses in our Milky 
Way and Andromeda models (2.6 x 10 12 Mq). This discrep- 
ancy suggests that the Milky Way and Andromeda do not 



contain the entire quantity of mass in the Local Group and 
are instead the most massive concentrations of mass within 
a larger all-encompassing medium. This point of view is 
very natural in a cosmological context where galaxies are 
not isolated islands, but rather mountain peaks within a 
vast continent of land. 

For these reasons our Local Group model supplements 
the Milky Way and Andromeda galaxy models with a diffuse 
and extended intragroup medium, as schematically depicted 
in Figure [T] Put within a cosmological framework this initial 
configuration may represent a point in time when the entire 
Local Group has decoupled from the general Hubble flow 
and is in the process of collapsing to become a virialized 
structure. 

For simplicity, we assume that the Local Group medium 
is initially a constant density cube of 1.5 Mpc on a side 
composed of both dark matter and gas. The total diffuse 
mass within this cube is set equal to the total mass of the 
two galaxy halos today, 2.6 x 10 12 Mq, yielding a net mass 
interior to the MW/M31 orbit that is consistent with the 
timing argument. O ur choice is mot ivated by other cosmo- 
logical simulations (|Gao et al.ll2004T l which indicate that a 
substantial fraction of the dark matter is likely to be diffuse 
and resides in between virialized halos (within the unviri- 
alized Local Group) at the initial time of our simulation 
~ 5 Gyr ago. 

We also postulate that the Local Group volume con- 
tains close to the cosmic mean value of baryons, namely 
16% (|Spergel et al.ll2003l ). Since the Milky Way and An- 
dromeda galaxy models only include baryons in the galactic 
disk and bulge components, they are far short of this value. 
We therefore set the Local Group medium to be 20% primor- 
dial gas, by mass, so that the entire region approaches the 
cosmic mean value. The gaseous component, like the dark 
matter, is initialized to have a constant density, and its tem- 
perature is fixed to be 3 x 10 5 K, consistent with estimates 
for the warm-hot intergalactic medium at comparable over- 
densities from co smological simulations (see, e.g., Fig. 6 in 
iDave et al.l l200lh and with observations of apparent RAM 
pressure stripping in the L ocal Group dwarf galaxy Pegasus 
|McConnachie et al.ll2007t ). 

We note that the initial gas temperature is far below 
the virialized temperature of the Local Group, however this 
assumption is consistent with the post-turn around, and 
pre-virialized initial conditions that we adopt. The constant 
density is also not in line with the assumed gr owth of th e ini- 
tial stellar disks from a virialized gaseous halo lMo et al.l (see, 
e.g.. Il998l . and references therein). Within a halo dynamical 
time, some of the initial intragroup gas is expected be ac- 
creted by each galaxy, and most of the remaining intragroup 
gas will end up being shock-heated during the collapse and 
relaxation of the Local Group. The dynamics of the merger 



3 The initial scale of the Local group in our simulation is an order 
of magnitude larger than the virial radius r2oo of the Milky Way 
or An dromeda galaxies. If one were to extend the iNavarro et all 
(1997) density profile of the envelope around each galaxy (with 
an asymptotic radial dependence of r~ 3 ) out to our initial Lo- 
cal Group scale, then one would roughly double the mass found 
within the virial radius of each halo. The added mass in that 
case would be similar to the amount we indeed assume for the 
intragroup medium. 
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Figure 1. Sketch of the initial configuration of our Local Group model, which consists of the Milky Way and Andromeda embedded in 
a diffuse, constant-density intragroup medium of equivalent mass. See the text in ij2]and Table 1 for more details. 



between the Milky Way and Andromeda is dominated by 
the distribution of dark matter and stellar disks, and so the 
basic results of our simulation, involving for example the 
timing of the merger, are not directly affected by the initial 
gas temperature and density distribution. 

Since the total mass of the diffuse medium is equiv- 
alent to the two galaxies, the two components are repre- 
sented with an equivalent number of particles; 1.3 million 
(1.04 million dark matter, and 260,000 gas). With these as- 
sumptions, the overdensity of our initial Local Group model 

iS <5lG = PLG//Ocrit ~ 10. 

2.4 Other Considerations 

Unfortunately, once the Local Group contains this extended 
mass distribution the relative motion of the Milky Way and 
Andromeda no longer becomes a trivial application of Ke- 
pler's laws. The diffuse mass will steadily extract orbital 
energy and angular momentum owing to dynamical friction, 
and the deep potential wells of the galaxies will slowly ac- 
crete diffuse matter. Both of these effects act to gradually 
extract energy from the binary system, hardening its orbit, 
and accelerating the merging process. 

Owing to the complicated nature of the orbit, and its 
deviation from simple two-body motion, we are left with a 
fair degree of ambiguity regarding the initial state of our 
Local Group model. Because the configuration envisaged by 
Figure [1] resembles the cosmological collapse of an overdense 
region of the Universe, a natural starting point is when this 
fluctuation decoupled from the general cosmological expan- 



sion, turned around, and began to collapse under its own self 
gravity. Like the original timing argument, we then presume 
that this scenario also occurred for the Milky Way and An- 
dromeda as well, i.e., we initialize their radial velocity to be 
zero as if they have just turned around and are now about 
to begin their gravitational collapse. 

While some insight about the value of the turn around 
radius can be gained from an application of the timing argu- 
ment, in practice we find that the position and velocity of the 
Milky Way and Andromeda quickly deviated from simple 
two-body motion owing to the diffuse intragroup medium. 
In general, the velocities quickly grew larger than the orig- 
inal orbit specified and the direction became highly radial. 
We therefore adopt a trial and error approach where we 
build an initial model, run it until the Milky Way and An- 
dromeda are separated by 780 kpc (as this is their current 
separation), and inspect the relative velocity between the 
two galaxies to assess the validity of the model. 

Employing this procedure results in a number of mod- 
els that are all honest representations of the Local Group 
at the present time. In practice, we find a trade off be- 
tween the initial separation and the initial orbital energy 
such that models which start with a larger separation re- 
quire a smaller initial eccentricity. In all cases, the apparent 
orbit of the two galaxies becomes (or was to begin with) 
nearly parabolic at the present time, consistent with re- 
cent estimates of the Milky Way- And romeda orbit (see, e.g., 
Ivan der Marel fc Guhathakurtall2007h . While the rest of this 
paper will primarily present the results of one particular 
model, we will also show that all models yield similar esti- 
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mates for the eventual merger between the Milky Way and 
Andromeda. We will argue in i]4.3.1l that this convergence 
results naturally from our assumed intragroup medium. 

The model we choose to focus upon begins with an 
initial separation of 1.3 Mpc, and initializes the Milky 
Way and Andromeda on an eccentric orbit e = 0.494, 
with a distance at perigalacticon of 450 kpc. With this or- 
bit the initial angular velocity is 65 km s~\ which could 
likely originate from tidal torques (|Gott fc Thuanl Il978l ; 
iRavchaudhurv fc Lvnden -Bcll 1989]). This particular model 
begins with the largest separation of all our models and 
therefore may be the best representation of the evolution 
of the Local Group since its decoupling from the univer- 
sal expansion. Since this model tracks the Local Group the 
farthest into the past, the intragroup medium also has a 
significant amount of time to react to the two galaxies and 
therefore is likely to be the most insensitive to its initial 
configuration. 



3 NUMERICAL METHODS 

To simulate the evolution of our Local Group and in partic- 
ular the interaction between the Milky Way and Andromeda 
use the pu blically availab le N-body/hydrodynamic code 
GADGET2 l|Springell l2005h . This version of the code em- 
ploys the "conservative-entrop y" formulation of Smoothed 
Particle Hydrodynamics (SPH, ISpringel fc Hernquistll2002l ) 
that conserves both en ergy and entropy (unlike earlier ver- 
sions of SPH; see e.g., lHernquisdll993br ). while improving 
shock-capturing. We assume that the gas is of primordial 
composition, and include the effects of radiative cooling. 

Star formation and its associated feedback are both 
included i n a m anner very similar to that described in 
ICox et ail (|2006l ). As is commonly assumed, stars are 
stochastically forme d at a rate determined by the SPH 



gas density (see, e.g . 
ISpringel et ail 120051 ; 



Katdll992l;ISpringel fc Hernquisdliooll ; 
Cox et alj|2006h with an efficiency set 



to match the ob served correlatio n between star formation 
and gas density (lKennicuttl ll998). 

Feedback from stellar winds and supernovae is treated 
in a very simplistic manner, namely the SPH particles that 
have sufficient density to form stars are fixed to have an ef- 
fective temperature of 10 5 K. This methodology is similar in 



back (see, e.i 


I.. Springel 2000 


; Springel & Hernquistl 12003; 


IStinson et al.ll2006l; Cox et al. 


200fjj), and is easy to imple- 



ment. Since the focus of this work is the large-scale evolution 
of the Local Group and the generic dynamics of the colli- 
sion between the Milky Way and Andromeda, the detailed 
treatment of the inter-stellar medium does not influence our 
primary conclusions. 

Numerical resolution is a significant consideration for 
any computational problem. For our purposes here, we re- 
quire sufficient resolution to reliably follow the interaction 
and merger of the Milky Way and Andromeda, while main- 
taining the ability to perform a number of simulations with 
the available computational resources. These considerations 
motivated the particle number choices outlined in fj2] Given 
the large number of components in these simulations (stellar 
disks, dark halos, and intragroup medium) and the desire to 
reduce two-body effects, we required all particles to have 



an identical mass of 2 x 10 7 Moand we employed a universal 
gravitational softening length of 150 pc. To test the sensi- 
tivity of the results to these parameter choices, we also ran 
a higher resolution version of one model with 30 times the 
baryonic disk mass resolution (and therefore number of par- 
ticles) and 2 times the resolution of the dark matter and 
intragroup medium. In this case we also decreased the grav- 
itational softening length of the baryons by a factor of 3, and 
increased that of the dark matter by a factor of 3. While this 
test yielded much better resolution of the stellar disks and in 
particular the tidal material, the general merger dynamics 
were identical to the low resolution version. 



4 THE COLLISION BETWEEN THE MILKY 
WAY AND ANDROMEDA 

4.1 Generic Features of the Merger 

In Figures through [5] we present the basic properties of 
the dynamical evolution of our Local Group, from 5 Gyr 
in the past and until 10 Gyr into the future, beyond the 
merger time between the Milky Way and Andromeda. Most 
of the features present in these figures are generic to bi- 
nary galaxy interactions, and hav e been described in great 
detail by prior studies (see, e.g ., iToomre fc Toomrej I l972 j; 
Barnes fc Hernq uist 1991, 1992; Mih os fc Hernquistl 11999 ; 



Cox et aJj|2006 j) . However, we will review some of the details 



that are particularly relevant to the Local Group, and sub- 
sequently highlight the unique status of our own Sun which 
will be a participant in this galaxy interaction. The future 
evolution of structures beyond the local group was simulated 
elsewhere (jNagamine fc Loebll200j . l2004l ; lBusha et allliool , 
l2005h . 

To begin, Figures[3]and[3]present the entire evolution of 
the Local Group from the point of view of a distant observer. 
These images begin at the start of our simulation, when 
the Milky Way and Andromeda are separated by 1.3 Mpc, 
and include the present state of the Local Group (labeled 
"Today") and the eventual merger of the Milky Way and 
Andromeda. As a guide to the eye, each panel includes the 
trajectory of both the Milky Way and Andromeda. 

Shown in Figure [5] is the evolution of the stellar com- 
ponent, which in our simulation only has contributions from 
the Milky Way and Andromeda as we ignore any structure 
smaller than the two largest galaxies in the Local Group. 

Figure |3] presents the projected gas distribution during 
the interaction, with panels shown at the same times as in 
Figure [5] Here the color-scale has been stretched to empha- 
size the abundant quantity of low-density gas that is spread 
throughout the local group. The initial condition of our Lo- 
cal Group model assumes a uniform distribution of warm 
gas, however the gas quickly responds to the non-uniform 
potential. In particular gas is accreted and shocked to form 
a hydrostatic halo of warm gas around the Milky Way and 
Andromeda galaxies. The gas distribution is also clearly af- 
fected by the interaction itself, as shocks develop once the 
galaxy halos begin to interpenetrate at close separation. 

While a detailed analysis of the diffuse gaseous intra- 
group medium is beyond the primary focus of our work, 
the basic properties of this component appear to be con- 
sistent with observations. In particular, our model pre- 
dicts that, at the present state, the intragroup medium is 
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Figure 2. Time sequence of the projected stellar density, shown in red-scale, for the merger of the Milky Way (MW) and Andromeda 
(M31). Andromeda is the larger of the two galaxies and begins the simulation in the upper-right. The Milky Way begins on the edge of 
the image in the lower-left. Each panel is 1.5 Mpc square, and the simulation time, in Gyr relative to today, appear on the top-left label 
of each panel. The trajectories of the Milky Way and Andromeda are depicted by the blue solid lines. 
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Figure 3. Same as Figure [2] only the projected gas density is shown. 
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Figure 4. Time sequence of the projected surface mass density of stars, during the final merger between Andromeda and the Milky 
Way. Panels have varying scales as specified by the label, in kpc, on the lower-right of each panel. The top-left panel is identical to the 
"Today" panel in Figure [2] The simulation time, with respect to today, is shown in the top-left of each panel. 
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predominantly warm 10 s — 10 6 K, has a fairly low den- 
sity, 10" 4 - 1(T 6 cm , and fills the entire volume that 
we simulate. We note this medium is also expected to 
extend to larger scales into what has been termed the 
"Warm-Hot Interga l actic Medium" (|Cen fc Ostrikerl 1 19991 ; 
iHellsten etail 1 19981 ; iDave et all 120011 ). While the presence 
of our intragroup m edium is consistent with all current data 
|Osone et al.l 12003 ). Owing to the difficulty in directly ob- 
serving gas in this state, most evidence for its existence 
comes from indirect means. For example, Chandra and 
FUSE observations of z ~ Oxygen and Neon absorption 
along numerous sight lines suggests th e presence of a local , 
volume-filling, diffuse, warm medium (|Nicastro et al.l |2002| . 
l2003l ; ISembacn et al.ll2003l ; ISavage et alj|2003h . although de- 
tailed analysis of the ionization states suggest that this is a 
complex multi-phase medium that our simulation does not 
have the resolution to model. 

While Figure[2]presented the dynamical evolution of the 
stellar mass on large scales, this vantage point makes it diffi- 
cult to distinguish the tell-tale signs of a galaxy merger. We 
therefore zoom into the central regions of the Local Group 
and specifically show the merger between the Milky Way 
and Andromeda in Figure [4] From this viewpoint the classic 
signatures of a galaxy interaction, such as tidal tails, plumes, 
and shells are clearly evident. 

The physical separation between the Milky Way and 
Andromeda is presented in Figure [5] The current state of 
the Local Group occurs ~ 5 Gyr after the start of our sim- 
ulation. For th e standard set of cosmological parameters 
|Tegmark| [2006). this implies that we initiate the simulation 
at a redshift of z ~ 0.5, around the time when the Sun was 
born in the Milky Way dislQ. Figure [5] also clearly shows the 
prediction that our models makes for the future collision be- 
tween the Milky Way and Andromeda. Their first close pas- 
sage will occur less than 2 Gyr from the present, and the cen- 
ters are fully coalesced in less than 5 Gyr. These time scales 
are c omparable to the lifetime of our Sun ((Sackman n et all 
1 19931 ) and admit the possibility that an observer in the Solar 
System will witness some (or all) of the galaxy collision. We 
will return to this possibility in H4.3.2I 

Finally, we present the relative velocity between the 
Milky Way and Andromeda in Figure [6] As in Figure [5] we 
clearly delineate the present state of the Local Group with 
a hatched region whose size corresponds to the errors esti- 
mated with the velocity measurements. We note that the 
velocities presented are relative to each galaxies center of 
mass, and do not correct for any motion relative to that. 

Now that the general features of our Local Group model 
have been outlined, we next explore the validity of our 
model, the merger dynamics, the star formation during the 
interaction, and the properties of the merger remnant. 



4.2 The Current State of the Local Group 




4 6 
Time (Gyr) 

Figure 5. The separation between the centers of Andromeda 
and the Milky Way during the course of their merger. The current 
separation of ~ 780 kpc is shown with a horizontal dashed line 
and occurs at T m 4.7 Gyr. 
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Figure 6. The relative velocity between the centers of An- 
dromeda and the Milky Way galaxies during the course of their 
merger. The present state of the Local Group occurs at T ss 
4.7 Gyr. 



Although our model for the Local Group presents its past, 
current, and future evolution, the only way to test its 



4 Note that extending the simulation to significantly earlier times 
is not adequate since stellar ages imply that the two galactic disks 
(and presumably their halos) have not been fully assembled at 
z > 2 jWvsdl2007t) . 



validity is by comparing it to the empirical data on its 
present-day state. As we focus primarily on the evolu- 
tion of the two largest galaxies in the Local Group, the 
Milky Way and Andromeda, it is only the relative sep- 
aration and motion of these two galaxies that can be 
compared to data. As mentioned in both iJTJ an d SI t ne 
separation and line-of-sight velocity of the Milky Way 
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and Andromeda are cur r ently measured to be 780 kpc 
(|McConnachie et alj|2005l ; iRibas et alj|2005l. and references 
therein), and -120 km s~ 1 l|Binnev fc Tremaind Il987n . re- 
spectively. The proper motion of Andromeda perpendicular 
to our line of sight is less well constrai ned, but current es- 
timates suggest t hat it is < 200 km s~ 1 |Peebles et al]l200ll ; 
lLoeb et alj|2005h . 

Figures [5] and [6] present the direct comparisons between 
the observational constraints and our model, and clearly 
demonstrate that it is viable. In particular, at a time of 
4.7 Gyr after the start of the simulation (close to the present 
cosmic time) the separation between the Milky Way and 
Andromeda is 780 kpc, while the relative radial velocity is 
135 km s _1 and the tangential velocity is 132 km s~ . While 
the tangential velocity is well within the limits currently fa- 
vored, the line-of-sight velocity is slightly larger than (but 
within 2.0 a of) the observations. Even though a number of 
our model assumptions can be manipulated to reduce these 
values, e.g., the initial separation and eccentricity of the 
fiducial orbit, we note that our model is a good fit to the 
velocity when the separation between the Milky Way and 
Andromeda is larger than 780 kpc. Given the observational 
uncertainties in these values we feel that there are likely to 
be a large number of models that can simultaneously fit all 
the data within its 2a error bars. 




-6 -4 -2 2 4 6 
Time Relative To Now (Gyr) 



Figure 7. The separation between the centers of Andromeda 
and the Milky Way during the course of their merger for a large 
ensemble of Local Group models. The current separation of ~ 
780 kpc is shown with a horizontal dashed line and all models are 
normalized to this particular time. 



4.2.1 An Ensemble of Models 



As mentioned in §231 the model described up to this point 
is one of twenty Local Group models that we have simu- 
lated. Since the primary conclusion of this paper, involving 
the timescale for the eventual merger between the Milky 
Way and Andromeda, may be influenced by any number of 
our model assumptions, we explicitly show the separation 
between the Milky Way and Andromeda and therefore the 
time of the merger, for all of our models in Figure [7] 

Figure [7] demonstrates several interesting features. 
Nearly all of the models provide a similar outcome for time 
when the two galaxies make their first passage (ensemble 
average and standard deviation are T = 2.8 ±0.5 Gyr). The 
same holds true for the final merger time (T = 5.4±0.4 Gyr). 
We also note that these average values are slightly larger 
than one provided by the model described up to this point. 
This highlights a general trend, namely the mergers that 
start with a larger separation and have to traverse a longer 
path through the intragroup medium, usually (but not al- 
ways) have a quicker merger dynamics. Regardless of the 
relatively small differences between merger times in these 
models, they are all completely coalesced by 6.2 Gyr from 
today. In the following section we will argue that this re- 
sult is a direct byproduct of our inclusion of an intragroup 
medium. 



4.3 Merger Dynamics 

4-3.1 Timescale 

One of the most intriguing characteristics of our Local 
Group model is the relatively quick timescale for the inter- 
action and merger between the Milky Way and Andromeda. 
As stated in the previous section, their first close passage 
will occur in less than 2 Gyr, and the final coalescence will 



occur in less than 5 Gyr. While the following section will 
specifically address the possible dynamics of the Sun during 
the merger, we show here that the cause for the short merger 
timescale is the intragroup medium. 

Our model for the Local Group assumes that the Milky 
Way and Andromeda are embedded in a diffuse medium of 
dark matter and gas. The medium itself (minus the evolved 
galaxies) is unvirialized and of low overdensity (S ~ 5) corre- 
sponding to a region of the Universe that has decoupled from 
the Hubble flow and started its evolution toward a collapsed 
virialized state. This medium exerts dynamical friction on 
the Milky Way and Andromeda and speeds up the merger 
dynamics by soaking orbital energy and angular momentum. 
This is shown explicitly in Figure [8] which compares the 
Milky Way - Andromeda separation in our fiducial model 
to an identical model without a diffuse medium. Without 
the intragroup medium, the merger timescale is nearly three 
times longer than with it. 

The effect of the intragroup medium should scale sim- 
ilarly to the standard (Chandrasekhar) formula for dynam- 
ical friction (|Binnev fc Tremaind Il9871 . Eq. 7-18), in which 
the deceleration of a massive object is proportional to the 
background matter density, dv/dt oc p. Therefore, the rate 
at which angular momentum is extracted from the orbit de- 
pends on the assumed intragroup medium density, a quan- 
tity which is poorly constrained observationally. Once the 
dark matter halos begin to interpenetrate, i.e., when the 
Milky Way - Andromeda separation is ~ 100 kpc, the 
merger completes relatively quickly because the dark matter 
halos dominate over the background density. 

This last point is particularly relevant to simulations of 
binary galaxy mergers, which typically omit any background 
overdensity. While this omission does not significantly al- 
ter the dynamical friction estimates once the halos overlap 
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Figure 8. The separation between Andromeda and the Milky 
Way during the course of their merger. Our fiducial model in 
red is compared to a case with no intragroup medium (dark 
matter+gas) in blue. Dynamical friction on the diffuse medium 
shortens the MW/M31 coalescence time by a factor of a few. The 
present time is 5 Gyr after the start of the simulation. 



and therefore dominate the background density, it may still 
change the distribution of orbits in high density environ- 
ments where the galaxies traverse through significant over- 
densities prior to their interaction event. In this sense, the 
merger timescales extracted from binary merger simulations 
may be considered an upper limit that is most applicable to 
galaxies in the field. 



4-3.2 The Fate of our Solar System 

An interesting consequence of the short timescale for the 
merger between the Milky Way and Andromeda is the pos- 
sibility that a human (or decedent thereof) observer will 
witness the interaction and merger. At the heart of this is- 
sue is the comparison between the lifetime of the Sun and 

the interaction timescale. 

Current evolutionary models (see, e.g. JSackmann et al.l 

1 19931 ) predict that the Sun will steadily increase its size 
and luminosity for the next 7 Gyr as it slowly consumes all 
available hydrogen and evolves towards a red giant phase. 
While this places a strong upper limit to the extent of life 
on Earth, it is likely that much smaller changes (< 50%) 
in the Sun's luminosity will significantly alter the Earth's 
atmosph ere and thus it s habitability within the next 1.1- 
3.5 Gyr llKastingl ll988r). lKorvcanskv et all ll200lD suggested 
that the onset of these effects could be delayed by increasing 
the orbital radius of the Earth through a sequence of inter- 
actions with bodies in the outer Solar System, and we can 
not rule out the possible colonization of habit able planets 



in ne arby stars, especially long-lived M-dwarf s (llldrv et al 



2007 ) whose lifetime may exceed 10 12 years (jAdams et al 



2005). In short, it is conceivable that life may exists for as 



little as 1.1 Gyr into the future or, if interstellar travel is 
possible, much longer. 

Regardless of the prospects for life in the future, we can 
attempt to predict what any potential observer at the solar 
Galactic circle might see by tracking candidate Suns in our 
simulation. In particular, we flag all stellar particles with a 
galactocentric orbital radius of 8±0.5 kpc, correspond ing to 
the observed value for the Sun l|Eisenhauer et al1l2003h . and 
subsequently follow the location of these particles forward 
in time. This procedure typically yields ~ 700 particles as 
"candidate Suns" (although it should be kept in mind that 
owing to our limited resolution, each of the simulated par- 
ticles is many orders of magnitude more massive than the 
Sun). The results are presented in Figure [9] and demonstrate 
some of the possible outcomes for the future location of our 
Sun. Given the uncertainties in our model parameters and 
the fact that the orbital period of the Sun around the Milky 
Way is much shorter than the merger timescale, it is not 
possible to forecast reliably the actual phase of the Galac- 
tic orbit of the Sun at the time of closest approach to An- 
dromeda. Therefore we regard all the stellar particles at the 
galactocentric radius of the Sun as equally probable of rep- 
resenting the Sun. We will first outline some of the general 
features of our fiducial model, before showing similar distri- 
butions from a subset of models to assess the reliability of 
these results. 

Figure [9] outlines the wide variety of potential locations 
for our candidate Suns during the future evolution of the 
Local Group. For example, the top-right panel in Figure [9] 
demonstrates the location of the candidate Suns after the 
first passage of Andromeda. At this point, the observer will 
most likely still be in the (now disturbed) disk of the Milky 
Way, but there exists a 12% chance that the Sun will be 
tidally ejected and take part in the tidal tail material that 
is > 20 kpc away from the Milky Way center. 

The probability that the candidate Sun will be located 
farther than 20 kpc from the Milky Way center steadily in- 
creases as the interaction progresses. At the second passage, 
the percentage of Suns that are farther than 20 kpc is 30%. 
This number increases to 48% after the second passage, and 
is 68% in the merger remnant. In fact, there is a 54% chance 
that the Sun will be at radii larger than 30 kpc in the merger 
remnant. However, we caution that this probability is de- 
rived at one point in time. Individual stars will generally 
spend much of their orbital time at large radii, even if their 
orbit is eccentric and they come much closer to the galactic 
center at other times. 

One unexpected possibility for the location of the fu- 
ture Sun is demonstrated in the lower-middle panel of Fig- 
ure [§] namely the Sun may actually become bound to An- 
dromeda instead of the Milky Way before the two galax- 
ies coalesce. Such a situation occurs when material becomes 
loosely bound after the first passage and is later captured 
by the gravitational potential of Andromeda during its sec- 
ond passage. While this outcome is unexpected and certainly 
exciting, only 2.7% of the candidate Suns became bound to 
Andromeda and so this outcome is relatively unlikely. 

In Figure [10] we attempt to assess whether the percent- 
ages just quoted for our fiducial model are characteristic 
of all our models and are therefore robust. Each panel in 
Figure [10] shows the distribution of galactocentric radii for 
candidate Suns during three periods of the interaction in 
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Figure 9. The possible location of the Sun during various stages of the merger between the Milky Way and Andromeda. The top panel 
shows all stellar particles in our simulation that have a present-day galactocentric radius of 8±0.1 kpc with a red cross and tracks their 
position into the future. The bottom panel presents a histogram of the radial distance from the center of the Milky Way. 



several representative Local Group models, including the 
fiducial one. The similarity between the models highlights 
the generic features of the fiducial model, although there 
is a large spread in the percentages quoted in the previ- 
ous paragraphs. The common distributions might have been 
expected given the similar orbital evolution demonstrated 
in Figure [7] and the tendency for mergers to preserve the 
hierarchy of the initial binding energies of the collisionless 
particles in the merger remnant. 

We note that several of the interactions generate unique 
tidal features which result in subtle differences in the distri- 
butions shown in Figure 1101 These differences are an out- 
come of different disk-spin orientations that were examined 
in the different models. One model in particular yields a rem- 
nant where some of the candidate Suns reside in a shell-like 



structure at ~ 50 — 60 kpc (see right-most panel in Fig- 
ure [10]). Only half of the models include the possibility for 
the Sun to be "stolen" by Andromeda, however one model 
yielded a probability that was nearly three times larger than 
in the fiducial case. 



4.4 Star Formation 

There is mounting evidence that galaxy interactions are the 
predominant mechanism for producing large bursts of star 
formation, s uch as in the ultra-luminous infrared galax- 
ies (U LIRGs, ISanders fc Mirabel|[l996l ; iBarnes fc Hernquistl 
1 1992T ) . During the interactions, gravitational torques extract 
angular momentum from the gas and funnel it to the center 
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Figure 10. The possible location of the Sun during the first passage (left panel), second passage (middle panel), and final merger (right 
panel) for several different models of the merger between the Milky Way and Andromeda. The histograms are difficult to distinguish, 
which highlights the generic probability for the Sun to reside at various galactocentric radii. The distribution of radii from the fiducial 
model, which is presented in Figure [9] is denoted by a hatched histogram. 




Time (Gyr) 



above the that which would occur if the Milky Way and An- 
dromeda has not participated in the merger at all. 

The weak starburst event triggered by the merger be- 
tween the Milky Way and Andromeda is a direct result of 
their present low gas content. Moreover, a large fraction 
(> 75%) of this gas will be consumed by quiescent star for- 
mation by the time the merger actually occurs. In short, 
both disks will be extremely gas-poor during the final coa- 
lescence and there will be no fuel for the starburst. 

While we have not explicitly tracked the black holes at 
the center of the Milky Way and Andromeda, it is interesting 
to speculate whether the merger will produce a luminous 
quasar, which many models ar gue is intricately lin ked to 
galaxy mergers and starbursts (Hopk ms et al.i r2006). Even 
though Figure [TJJ demonstrates that there is not enough 
gas to fuel a powerful starburst, this gas content is clearly 
sufficient to ignite a luminous quasar if ~ 1% of it is accreted 
by the black hole. While the current work cannot address 
this possibility in detail, our model provides a framework to 
study the formation of quasars in the future. 



Figure 11. The cumulative star- formation rate during the 
merger of the Milky Way and Andromeda compared to the star 
formation for models of the Milky Way and Andromeda evolved 
in isolation. 



of the merger galaxy where it participates in a centrally- 
concentrated starburst. Because the future of the Local 
Group entails a major galaxy merger, it is natural to exam- 
ine whether the merger of the Milky Way with Andromeda 
will become a ULIRG. 

This question is addressed explicitly in Figure [11] which 
shows the star-formation rate during the entire evolution of 
the Local Group, including the merger between the Milky 
Way and Andromeda. Throughout the entire evolution of 
the Local Group, the star-formation rate steadily decreases 
with time, and hence we conclude that the merger Local 
Group will not become a ULIRG in the future. In fact, the fi- 
nal coalescence yields star formation that is barely enhanced 



4.5 The Merger Remnant 

Galaxy mergers have become an area of intense study ow- 
ing to their proposed role in shaping g alaxy morpho logy. In 
particular, the "merger hypothesis" I Toomrel Il977l ) posits 
that the interaction and merger of two spiral galaxies leaves 
behind a remnant that is morphologically and kinematically 
similar to an elliptical galaxy. Since the future evolution of 
the Local Group contains such an event, it is natural to 
examine whether the local group will eventually consist of 
a single elliptical galaxy, and if so, how do the properties 
of this galaxy differ from present-day ellipticals which were 
formed at earlier epochs. 

Note that the morphology of the Milky Way and An- 
dromeda merger remnant, which we abbreviate hereafter as 
Milkomeda, was presented in Figure [4] This figure confirms 
the notion that galaxy mergers leave behind remnants that 
are spheroidal in shape, contain stars with a wide range of 
orbits and a large velocity dispersion, and possess very mod- 
est net rotation. 
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Figure 12. Properties of the Milky Way and Andromeda merger remnant, which we abbreviate as Milkomeda. Panel (a) shows the 
spherically averaged distribution of dark matter, stars, and gas. The mass density is expressed in units of the present-day critical density, 
Pcrit = 10~ 2<J g cm -3 . Panel (b) shows the projected distribution of stellar mass, including the contribution from stars formed during 
the course of the simulation which are termed "New Stars," plotted against i? 1 / 4 , where R is the projected radius. Panel (c) presents the 
kinematics of Milkomeda. For 195 projections, the ellipticity of the half-mass isophote e, the rotation velocity (V) along the major axis 
V ma j and the central velocity dispersion er, are plotted as a colored histogram, with darker shades representing a higher concentration of 
projected properties. The solid line represents what is expected from an oblate isotropic rotation and overplotted is data extracted from 
the literature. Panel (d) shows the velocity anisotropy, defined as the ratio V/cr normalized by the value expected for an oblate isotropic 
rotator, versus the isophotal shape parameter 04/0, where positive values are disky isophotes and negative values are boxy. Overplotted 
in black are current data (Ralf Bender, private communication). Panel (e) shows th e X-ray surface brig htness profile, and panel (f) shows 
Milkomeda on the near-IR fundamental plane along with the relation observed bv lPahre et al 
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In Figure 1121 we present a series of plots that further 
support the assertion that Milkomeda resembles an elliptical 
galaxy and also serve to better quantify its properties. In 
panel (a), we present the spherically averaged mass profile, 
decomposed by the three components; dark matter, stellar 
and gaseous mass. At radii greater than ~ 2 — 3 kpc, the dark 
matter dominates the mass density, and at ra dii greater than 
~ 20 kpc, the profil e is well-fit by th e NFW (jNavarro et al.l 
119961 ) or Hernquist (|Hernquistlll990l ) profile. 

The projected mass distribution shown in panel (b) of 
Figure [T^] presents the first direct evidence that Milkomeda 
resembles an elliptical galaxy. Specifically, the stellar sur- 
face density is close to a pure R 1//4 distribution, with the 
exception of the inner ~ 400 pc (~ 2.5 x the gravitational 
softening length), where the surface density flattens to a 
nearly constant density core. 

We have also quantified the kinematics and isopho- 
tal shape of t he M ilkomeda galaxy following the methods 
of ICox et all <|2006h . Panel (c) in Figure IT21 presents the 
anisotropy diagram, a measure of the half-mass isophote 
ellipticity versus the maximum rotation along its major axis 
divided by the central velocity dispersion. The shaded re- 
gion represents the distribution of values for Milkomeda if 
viewed from 195 directions that uniformly sa mple the unit 
spher e (with angles selected using HEALPIX, iGorski et al.1 
120051 ). Also plotted in this figure is the relation expected 
from an oblate isotropic rotator as a solid line. Owing to the 
high concentration of projections that closely track the solid 
line, this analysis demonstrates that Milkomeda is nearly an 
oblate isotropic rotator. In addition, the deviations from 
a perfect ellipse are also quantified and are presented in 
panel (d) in Figure [T2"l Depending on the viewing direction, 
Milkomeda may appear to be either disky or boxy, with a 
slightly larger fraction of views producing disky isophotes. 

We have also analyzed the properties of the hot gas 
in and around Milkomeda. The evolution of this component 
was presented in Figure [3] and clearly demonstrates the for- 
mation of an extended gaseous halo, primarily accreted from 
the large reservoir of the intragroup medium. Although the 
gas temperature was originally 3 x 10 5 K, it has been shock- 
heated to the virial temperature of ~ 3 x 10 6 K in Milkomeda 
and has a gradient to cooler temperatures at large radii. This 
hot gas leads to an X-ray surface brightness profile shown 
in panel (e) of Figure [12] and a total X-ray luminosity of 
~ 10 41 ergs s _1 , which is consistent with present-day ellipti- 
cal galaxies of equivalent B-band luminosity (~ 3x 10 10 Lq). 

In general, Milkomeda resembles the remnants of gas- 
rich major mergers, which in-turn resemble the general pop- 
ulation of low- and moderate-lum inosity elliptical galaxies 
(|Naab et alj200d : ICox et al.|[2006h . However, there are some 
systematic differences that likely arise because of the much 
smaller gas content of the Milky Way and Andromeda when 
they merge. In particular, the inner regions of Milkomeda 
have a much lower stellar density than present-day ellipti- 
cals, which are often observed to have ste e p spikes of newly 
form ed stars l|Rothberg fc Joseph! |2004| ; iKormendv et al.l 
120071 ). These central excesses arise from high gas concentra- 
tions th at fuel nuclear starbursts during the merger event 
(see e .g., lMihos fc Herncmistll 19941 ; ISpringel|[2000l ; ICox et all 
2006). This process does not occur during the formation of 
Milkomeda (see £|4.4f> . for which the stars formed during the 



merger simulation have an identical profile to the entire stel- 
lar population (see panel (b) in Figure [12)) . 

The diffuse nature of Milkomeda is also evident if it 
is projected onto the near-IR fundamental plane as shown 
in panel (f) of Figure 1121 where it lies above the observed 
relation. The half-mass radius is 4.9 kpc, which is larger 
than the mean relation found in the Sloan D igital Sky Sur- 
vey l|Shen et al.ll2003l ; iDesroches et al.l l2007l for galaxies of 
equivalent stellar mass (1.3 x 1O 12 M0) or r-band absolute 
magnitude (-21.2). These comparisons gives credence to the 
claim that present-day ellipticals can not have formed fro m 
the merger of present day spirals ( Na ab fc Ostrik cr 2007). 



5 CONCLUSIONS 

In this paper we have used an N-body/hydrodynamic sim- 
ulation to track the evolution of the Local Group, focusing 
primarily upon the two most massive galaxies: the Milky 
Way and Andromeda. In contrast to most prior work, which 
typically employed models for the Local Group to infer its 
total mass or the proper motion of its constituents, we sim- 
ulated the large-scale dynamics of all matter in the Lo- 
cal Group including dynamical friction on the intragroup 
medium, leading to the eventual merger of the Milky Way 
and Andromeda. 

Owing to the diffuse intragroup medium that was as- 
sumed to pervade the Local Group with a total mass compa- 
rable to that the galaxies, we have found that the interaction 
and merger between the Milky Way and Andromeda will oc- 
cur in less t han 5 Gyr, a timescale comparable to the lifetime 
of the Sun (|Sackmann et al]|l993T ) . This Local Group model 
therefore admits the possibility that future astronomers in 
the Solar System will witness parts of, or the entire interac- 
tion and merger of the Milky Way and Andromeda. 

With this in mind we have calculated the probable lo- 
cation of our Solar System during specific points of the fu- 
ture interaction and find several interesting outcomes. First, 
there is a chance that the Sun will be ejected along with 
other tidal material into a long tidal tail following the first 
passage of Andromeda. Second, as a result of the disrup- 
tive effects of each close tidal passage, there is an increasing 
chance that the Sun will inhabit extended tidal features as 
the interaction proceeds. Moreover, there is a small chance 
that the Sun will be more tightly bound to Andromeda at 
some point during the merger. In such a case, Andromeda 
will capture the Sun and future astronomers in the solar sys- 
tem might see the Milky Way as an external galaxy in the 
night skjQ. 

While this paper highlights the possible outcomes of the 
future interaction between the Milky Way and Andromeda, 
we emphasize that our model is likely to be one of many 
plausible models within an ensemble of possibilities that 



5 Aside from the changed night sky during the interaction be- 
tween the Milky Way and Andromeda, future observers might 
witness enhanced comet showers due to the increased flux of stars 
passing by the solar system and perturbing the Oort cloud. Ad- 
ditional observables would include a slightly enhanced star for- 
mation in the two galaxies and the production of hot gas by the 
shocks surrounding the interaction region. 
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span the uncertain value of the transverse velocity of An- 
dromeda and the density (and other properties) of the in- 
tragroup medium. We have performed 20 additional runs in 
order to test the sensitivity of our results to various assump- 
tions of our model - mainly involving the initial orbit of the 
Milky Way and Andromeda. These runs yield similar esti- 
mates for the merger timescale as well as for the possible lo- 
cations of the Sun in the future, provided that the intragroup 
medium is indeed similar to our fiducial case. While this 
gives us some confidence that our results are robust, an even 
larger suite of models, that spans a much wider set of model 
assumptions, will provide better statistics on these results. 
In addition, employing higher resolution simulations with 
increased complexity, may shed light on the nat ure of the in- 
tragr oup medium, the soft X -ray background (lOsone et al.l 
2002), galactic substructure |Willman et al.l [20051) . the or i- 
gin of the Magellanic Clouds and Stream (|Besla et al.ll2007l). 
and the future evolution of globular clusters (jForbes et al.l 
l2000h . 

Finally, we note that the simulated views from the dis- 
tribution of locations for the candidate Suns in the merger 
remnant (see Fig. [9]), which we have termed Milkomeda, rep- 
resent the only views available for a future local astronomer. 
Extragalactic astronomy will come to an end within 100 bil- 
lion years if the cosmological constant will not evolve with 
time. Owing to the accelerated expansion caused by a steady 
cosmological constant, all galaxies not bound to the Local 
Group will eventually rece de away fro m the Local Group 
and exit our event horizon |Loeb| [2002). At that point, the 
merger product of the Milky Way and Andromeda (with its 
bound satellites) will con stitute the entire visible Universe 
l|Nagamine fc Loebl [20031 ). 
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